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Abstract. We are introducing a model-free control and a control with a restricted 
model for finite-dimensional complex systems. This control design may be viewed 
as a contribution to "inteUigent" PID controllers, the tuning of which becomes 
quite straightforward, even with highly nonlinear and/or time-varying systems. 
Our main tool is a newly developed numerical differentiation. Differential algebra 
provides the theoretical framework. Our approach is validated by several numerical 
experiments. 



1. INTRODUCTION 

Writing down simple and reliable differential 
equations for describing a concrete plant is al- 
most always a daunting task. How to take into 
account, for instance, frictions, heat effects, ageing 
processes, characteristics dispersions due to mass 
production, . . . ? Those severe difficulties explain 
to a large extent why the industrial world is 
not willing to employ most techniques stemming 
from "modern" control theory, which are too often 
based on a "precise" mathematical modeling, in 
spite of considerable advances during the last fifty 
years. We try here ^ to overcome this unfortunate 



^ This communication is a slightly modified and updated 
version of (Fliess & Join [2008a]), which is written in 
French. Model-free control and control with a restricted 
model, which might be useful for hybrid systems (Bourdais, 
Fliess, Join & Perruquetti [2007]), have already been 
applied in several concrete case-studies in various domains 
(Choi, d'Andrea-Novel, Fliess & Mounier 12009], Gedouin, 
Join, Delaleau, Bourgeot, Chirani & Calloch 12008], Join, 



situation thanks to recent fast estimation meth- 
ods. ^ Two cases are examined: 

(1) Model-free control is based on an elementary 
continuously updated local modeling via the 
unique knowledge of the input-output behav- 
ior. It should not be confused with the usual 
"black box" identification (see, e.g., (Ker- 
schen, Worden, Vakakis & Golinval [2006], 
Sjoberg, Zhang, Ljung, Benveniste, Delyon, 
Glorennec, Hjalmarsson & .luditsky [1995])), 
where one is looking for a model which is 
valid within an operating range which should 
be as large as possible. ^ Let us summarize 
our approach in the monovariable case. The 



Masse & Fliess 12008], Villagra, d'Andrea-Novel, Fliess & 
Mounier l2008a,b]). Other applications on an industrial 
level are being developed. 

^ See, e.g., (Fliess, Join & Sira-Ramirez 12008], Mboup, 
Join & Fliess 12009]) and the references therein. 
^ This is why we use the terminology "model-free" and not 
"black box". 



input-output behavior of the system is as- 
siuiiod to "approximativcly" governed within 
its operating range by an unknown finite- 
dimensional ordinary differential equation, 
which is not necessarily linear, 

i;(2/,y,...,y(«),n,w,...,n(^)) = (1) 

We replace Eq. (1) by the following "pheno- 
menological" model, which is only valid dur- 
ing a very short time interval, 

y(^) =F + au (2) 

The derivation order v, which is in general 
equal to 1 or 2, and the constant parameter a 
are chosen by the practitioner. It implies that 
V is not necessarily equal to the derivation 
order a of y in Eq. (1). The numerical value 
of F at any time instant is deduced from 
those of u and y'^"^ , thanks to our numerical 
differentiators (Fliess, Join & Sira- Ramirez 
[2008], Mboup, Join & Fliess [2009]). The de- 
sired behavior is obtained by implementing, 
if, for instance, v = 2, the intelligent PID 
controller'^ {i-PID) 

u = -- + ^+ Kpe + Ki [ e + Kdb (3) 
a a J 

where 

• y* is the output reference trajectory, 
which is determined via the rules of 
fiatness-based control (see, e.g., Fliess, 
Levine, Martin & Rouchon [1995], Rotella 
& Zambettakis [2007], Sira- Ramirez & 
Agrawal [2004]); 

• e = y — y* is the tracking error; 

• Kp, Kj, Ku are the usual tuning gains. 
(2) Assume now that a restricted or partial 

model of the plant is quite well known and 
is defined by Eq. (1) for instance. The plant 
is then governed by the restricted, or incom- 
plete, modeling, ^ 

E{y, ij,..., y'-''\u, u,..., u*^)) +G = 

where G stands for all the unknown parts. ^ 
If the known system, which corresponds to 
= 0, is fiat, we also easily derive an 
intelligent controller, or i-controller, which 
gets rid of the unknown effects. 

Differential algebra is briefiy reviewed in Sect. 2 
in order to 

• derive the input-output differential equa- 
tions. 



* This terminology, but with other meanings, is not new in 
the literature (see, e.g., Astrom, Persson & Hang [1992]). 
^ Since we are not employing the terminology "black box", 
we are also here not employing the terminology "grey box". 
® Any mathematical modeling in physics as well as in 
engineering is incomplete. Here the term G does not need 
to be "small" as in the classic approaches. 



• define minimum and non-minimum phase 
systems. 

Numerical differentiation of noisy signals is exam- 
ined in Sect. 3. Sect. 4 states the basic princi- 
ples of our model-free control. Several numerical 
experiments^ arc reported in Sect. 5. We deal 
as well with linear® and nonlinear systems and 
with monovariable and multivariable systems. An 
anti-windup strategy is sketched in Sect. 5.4.2. 
Sect. 6 studies the control with a partially known 
medehng. One of the two examples deals with 
a non-minimum phase system. ^ The numerical 
simulations of Sect. 5.1 and Sect. 6.2 show the that 
our controllers behave much better than classic 
PIDs. Several concluding remarks are discussed 
in Sect. 7. 

Remark 1. Only Sect. 2 is written in an abstract 
algebraic language. In order to understand the 
sequel and, in particular, the basic principles of 
our control strategy, it is only required to admit 
the input-output representations (1) and (4), as 
well as the foundations of fiatness-based control 
(Fliess, Levine, Martin & Rouchon [1995], Rotella 
& Zambettakis [2007], Sira- Ramirez & Agrawal 
[2004]). 

2. NONLINEAR SYSTEMS 
2. 1 Differential fields 

All the fields considered here are commutative and 
have characteristic 0. A differential field^^ ^ is a 
field which is equipped with a derivation ^, i.e., 
a mapping M. such that, y a,b £ A, 

• ^{a + b) =d + b. 



^ Those computer simulations would of course be impos- 
sible without precise mathematical models, which are a 
priori known. 

® The remark 9 underlines the following crucial fact: the 
usual mathematical criteria of robust control are becoming 
pointless in this new setting. 

^ Non-minimum phase systems are today beyond our 
reach in the case of model-free control. This is certainly 
the most important theoretical question which is left open 
here. 

^"We are perfectly aware that such a comparison might be 
objected. One could always argue that an existing control 
synthesis in the huge literature devoted to PIDs since 
Ziegler & Nichols [1942] (see, e.g., Astrom & Ilagglund 
]2006], Astrom & Murray ]2008], Besangon-Voda & Gentil 
[1999], Dattaet, Ho & Bhattacharyya 12000], Dindeleux 
11981], Franklin, Powell & Emami-Naeini 12002], John- 
son & Moradi 12005], Lequesne [2006[, O'Dwyer [2006[, 
Rotella & Zambettakis [2008], Shinskey [1996], Visioli 
[2006], Wang, Ye, Cai & Hang [2008], Yu [1999]) has been 
ignored or poorly understood. Only time and the work of 
many practitioners will be able to confirm our viewpoint. 
"See, e.g., (Chambert-Loir [2005], Kolchin [1973[) for 
more details and, in particular, (Chambert-Loir [2005]) for 
basic properties of usual fields, i.e., non-differential ones. 



• ^(ab) = db + ab. 

A constant c G ^ is an element such that c = 0. 
The set of all constant elements is the subfield of 

constants. 

A differential field extension is defined by two 
differential fields K, £ such that: 

• the derivation of ^ is the restriction to ^ of 
the derivation of £. 

Write ^{S), S c £, the differential subfield of £ 
generated by A and S. Assume that £,/A is finitely 
generated, i.e., £ = A{S), where S est finite. An 
element ^ G £ is said to be differentially algebraic 
over R if, and only if, it satisfies an algebraic 
differential equation . . . ,^'"^) = 0, where P 
is a polynomial function over .ft in n+ 1 variables. 
The extension £/.ft is said to be differentially 
algebraic if, and only if, any element of £ is 
differentially algebraic over R. The next result is 
important: 

£,/M is differentially algebraic if, and only if, its 
transcendence degree is finite. 

An element of £, which is non-differentially al- 
gebraic over is said to differentially transcen- 
dental over R. An extension £/.ft, which is non- 
differentially algebraic, is said to be differentially 
transcendental. A set {^t € £ | /. S /} is said to be 
differentially algebraically independent over R if, 
and only if, there does not exists any non-trivial 
differential relation over R: 

Q(...,e^\...) = 0, where Q is a polynomial 
function over R, implies Q = 0. 
Two such sets, which are maximal with respect 
to set inclusion, have the same cardinality, i.e., 
they have the same number of elements: this is the 
differential transcendence degree of the extension 
£/.ft. Such a set is a differential transcendence 
basis. It should be clear that £,/R is differentially 
algebraic if, and only if, its differential transcen- 
dence degree is 0. 

2.2 Nonlinear systems 

2.2.1. General definitions Let k he & given dif- 
ferential ground field. A system^"^ is a finitely 
generated differentially transcendental extension 
IC/k. Let m be its differential transcendence de- 
gree. A set of (independent) control variables u = 
{ui,...,Um) is a differential transcendence ba- 
sis of /C/fc. The extension IC/k{u) is therefore 
differentially algebraic. A set of output variables 
y = {yi, . . . , yp) is a subset of /C. 



12 See also (Delaleau [2002, 2008], Fliess, Join & Sira- 
Ramirez [2008], Fliess, Levine, Martin & Rouchon ]1995]) 
which provide more references on the use of differential 
algebra in control theory. 



Let n be the transcendence degree of K/k{u) and 
let X = {xi, . . . ,Xn) be a transcendence basis 
of this extension. It yields the generalized state 
representation: 

«,...,«(")) =0 
BK{yK,x,u,...,u'^^^) = 

where A^, l = l,...,n, B^,, k — l,...,p, are 
polynomial functions over k. 

The following input-output representation is a 
consequence from the fact that yi, . . . ,yp are 
differentially algebraic over k{u): 

<I>,(y,...,y(^^),«,...,w(*^))-0 (4) 

where j = 1,. . . ,p, is a polynomial function 
over k. 

2.2.2. Input- output invertibility 

• The system is said to be left invertible if, 
and only if, the extension k{u,y)/k{y) 
is differentially algebraic. It means that one 
can compute the input variables from the 
output variables via differential equations. 
Then m< p. 

• It is said to be right invertible if, and only 
if, the differential transcendence degree of 
k{y)/k is equal to p. This is equivalent saying 
that the output variables are differentially 
algebraically independent over k. Then p < 
m. 

The system is said to be square if, and only if, 
m = p. Then left and right invertibilities coincide. 
If those properties hold true, the system is said to 
be invertible. 

2.2.3. Minimum and non-minimum phase sys- 
tems The ground field k is now the field M of real 
numbers. Assume that our system is left invert- 
ible. The stable or unstable behavior of Eq. (4), 
when considered as a system of differential equa- 
tions in the unknowns u {y is given) yields the 
definition of minimum, or non-minimum, phase 
systems (compare with Isidori [1999]). 



3. NUMERICAL DIFFERENTIATION 

The interested reader will find more details and 
references in (FHess, Join & Sira- Ramirez [2008]). 
We refer to (Mboup, Join & Fliess [2009]) for 
crucial developments which play an important role 
in practical implementions. 



-"^^This extension k{u,y)/k{y) is called the residual dy- 
namics, or the zero dynamics (compare with Isidori ]1999]). 



3. 1 General principles 



4-1 Local modeling 



3.1.1. Polynomial signals Consider the polyno- 
mial time function of degree N 



N 



where t>0. Its operational, or Laplace, transform 
(see, e.g., Yosida [1984]) is 



(5) 



v=0 



Introduce j^, which is sometimes called the alge- 
braic derivation. Multiply both sides of Eq. (5) by 
^s^+i, a = 0,l,...,iV. The quantities x<^''\0), 
V = 0,1,..., A^, which satisfy a triangular sys- 
tem of linear equations, with non-zero diagonal 
elements. 



N 



^ y 

\v=0 



(0)s 



JV- 



(6) 



are said to be linearly identifiable (Fliess & Sira- 
Ramirez p003]). One gets rid of the time deriva- 
tives si^^^, fj, = 1,...,N, < i < TV, by 
multiplying both sides of Eq. (6) by , N > N. 

Remark 2. The correspondence between and 
the product by (— t)" (sec, e.g., Yosida [1984]) 
permits to go back to the time domain. 



(1) Assume that the system is left invertible. If 
there are more output variables than input 
variable, i.e., p ^ m, pick up m output 
variables, say the first m ones, in order to 
get an invertible square system. Eq. (2) may 
be extended by writing 

y("i) = Fi+ ai,iUi H h ai^rnUm 

(7) 

(n ) 

Up " = Fm-\- am,lUi H 1- ap^mUm 

where 

• nj > 1, i = l,...,p, and, most often, 
Uj = 1 or 2; 

• aj,i e K, i = 1, . . . ,m, j — 1,. . . ,p, are 
non-physical constant parameters, which 
are chosen by the pratitioner such that 
aj^iUi and Fj are of the same magnitude. 

(2) In order to avoid any algebraic loop, the 
numerical value of 

= y) ~ aj,mUm 

is given thanks to the time sampHng 

m 



1=1 



^iUi{K - 1) 



where [•('«)]e stands for the estimate at the 
time instant k. 
(3) The determination of the reference trajecto- 
ries for the output variables yj is achieved in 
the same way as in flatness-based control. 



3.1.2. Analytic signals A time signal is said to 
be analytic if, and only if, its Taylor expansion is 
convergent. Truncating this expansion permits to 
apply the previous calculations. 



3.2 Noises 

Noises are viewed here as quick fluctuations 
around 0. They are therefore attenuated by low- 
pass filters, like iterated integrals with respect to 
time. 



4. MODEL-FREE CONTROL: GENERAL 
PRINCIPLES 

It is impossible of course to give here a complete 
description which would trivialize for practitioners 
the implementation of our control design. We hope 
that numerous concrete applications will make 
such an endeavor feasible in a near future. 



'^^See (Fliess [2006]) for a precise mathematical theory, 
which is based on nonstandard analysis. 



Remark 3. It should also be pointed out that, in 
order to avoid algebraic loops, it is necessary that 
in Eq. (4) 

^0, j = l,...,p 



It yields 



rij < Nj 



94>, 



Numerical instabilities might appear when — 
is closed to 0. 

Remark 4- Our control design lead with non- 
minimum phase systems to divergent numerical 
values for the control variables Uj and therefore 
to the inapphcability of our techniques. 

Remark 5. Remember that the Equations (1) and 
(4) are unknown. Verifying therefore the proper- 
ties discussed in the two previous remarks may 
only be achieved experimentally within the plant 
operating range. 



^^This kind of difficulties has not yet been encountered 
whether in the numerical simulations presented below nor 
in the quite concrete applications which were studied until 
now. 



4-2 Controllers 

Let us restrict ourselves for the sake of notations 
simplicity to monovariable systems. If i/ = 2 in 
Eq. (2), the intelligent PID controller has already 
been defined by Eq. (3). If = 1 in Eq. (2), replace 
Eq. (3) by the intelligent PI controller, or i-PI, 



F y* 

u= h — + Kpe + Ki 

a a 



(8) 



Remark 6. Until now we were never obliged to 
chose V ^ 2 m Eq. (2). The previous controllers 
(3) and (8) might then be easily extended to to 
the generalized proportional integral controllers, 
or GPIs, of (Fhess, Marquez, Delaleau & Sira- 
Rami'rez [2002]). 

Remark 7. In order to improve the performances 
it might be judicious to replace in Eq. (3) or in 
Eq. (8) the unique integral term Kj J e by a finite 
sum of iterated integrals 



where 

• J ... J e stands for the iterated integral of 

order A, 

• the Ki^ , A = 1, . . . , A, are gains. 

We get if Ki^ ^ an intelligent PI^D or PI^ 
controller. Note also that setting A = is a 
mathematical possibility. The lack of any integral 
term is nevertheless not recommended from a 
practical viewpoint. 

Let us briefly compare our intelligent PID con- 
trollers to classic PID controllers: 

• We do not need any identification procedure 
since the whole structural information is con- 
tained in the term F of Eq. (2), which is 
eliminated thanks to Eq. (3). 

• The reference trajectories, which are chosen 
thanks to flatness-based methods, is much 
more flexible than the trajectories which are 
usually utihzed in the industry. Overshoots 
and undershoots are therefore avoided to a 
large extent. 



5. SOME EXAMPLES OF MODEL-FREE 
CONTROL 

A zero-mean Gaussian white noise of variance 0.01 
is added to all the computer simulations in order 
to test the robustness property of our control 
design. We utilize a standard low-pass fllter with 



classic controllers and the principles of Sect. 3.2 
with our intelligent controllers. 



5.1 A stable monovariable linear system 



The transfer function 



(^ + 2)^ 

(S + 1)3 

deflnes a stable monovariable linear system. 



(9) 



5.1.1. A classic PID controller We apply the 
well known method due to Broi'da (see, e.g., 
Dindeleux [1981]) by approximating the system 
(9) via the following delay system 



{Ts+1) 

K = 'i,T = 2.018, T = 0.2424 are obtained thanks 
to graphical techniques. The gain of the PID 
controller are then deduced (Dindeleux [1981]): 

^ 100(0. 4r+T) _ -, n-in-i _ 1 _ 

0.7754, Kd = 2^ = 0.1766. 



5.1.2. i-PI. We are employing y — F + u and 
the i-PI controller 



u = -[Fl + y* + Pl{e) 



where 



• [F]e = [y]e - U, 

• is a reference trajectory, 

• e = y-y*, 

• PI(e) is an usual PI controfler. 

5.1.3. Numerical simulations Fig. 1 shows that 
the i-PI controller behaves only slightly better 
than the classic PID controller. When taking into 
account on the other hand the ageing process and 
some fault accommodation there is a dramatic 
change of situation: 

• Fig. 2 indicates a clear cut superiority of our 
i-PI controller if the ageing process corre- 
sponds to a shift of the pole from 1 to 1.5, 
and if the previous graphical identification is 
not repeated. 

• The same conclusion holds, as seen Fig. 3, if 
there is a 50% power loss of the control. 

Remark 8. This example shows that it might use- 
less to introduce delay systems of the type 



T{8)e 



-Ls 



reiR(s), L>o 



(10) 



^^The extension to the multivariable case is immediate. 
See Sect. 5.3. 



for tuning classic PID controllers, as often done 
today in spite of the quite involved identification 
procedure. It might be reminded that 



• the structure and the control of systems 
of type (10) have been studied in (FUess, 
Marquez & Mounier [2002]), 

• their identification with techniques stemming 
also from (FHess & Sira- Ramirez [2003]) has 
been studied in (Belkoura, Richard & FHess 
[2009], Ollivier, Moutaouakil & Sadik [2007], 
Rudolph & Woittennek [2007]). 

Remark 9. This example demonstrates also that 
the usual mathematical criteria for robust control 
become to a large irrelevant. Let us however 
point out that our control leads always to a pure 
integrator of order 1 or 2, for which the classic 
frequency techniques (see, e.g., Astrom & Murray 
[2008], Franklin, Powell & Emami-Naeini [2002], 
Rotella & Zambettakis [2008]) might still be of 
some interest. 

B.ernark 10. As also shown by this example some 
fault accommodation may also be achived without 
having recourse to a general theory of diagnosis. 

5.2 A monovariable linear system with a large 
spectrum 

With the system defined by the transfer function 



(s + l)(s + 0.1)(s + 0.01)(s - 0.05)(s - 0.5)(s - 5) 

We utilize y = F+u. A i-PI controller provides the 
stabilization around a reference trajectory. Fig. 4 
exhibits an excellent tracking. 

5.3 A multivariable linear system 
Introduce the transfer matrix 

f (s + 0.01)(s + 0.1)(s-l)^, ° J 

\ (s+0.003)(s-0.03)(s + 0.3)(s + 3) (,s + 0.004) (s+0.04) (s-0.4) (s + 4) / 

For the corresponding system we utilize after a 
few attempts Eq. (7) with the following decoupled 
form 

in = Fi+ IOmi 1)2 =F2 + 1{)U2 

The stabilization around a reference trajectory 
(2/1,2/2) is ensured by the multivariable i-PID 
controller 

"1 = 1^ (^1 Kpiei + Kii / ei + Koiei) 

U2 = Jo (j/2 -F2+ Kp2e2 + K12 / 62 + K 02(^2) 

where 

• ei = yt-yi, 62 = 1/2-2/2; 

• Kpi = 1, Kn = Kdi = 0, Kp2 = K12 = 50, 

Kd2 = 10. 

The performances displayed on Fig. 5 and 6 are 
excellent. Fig. 6-(b) shows the result if we would 
set Fi = i^2 = 0: it should be compared with Fig. 
6-(a). 



Remark 11. Model reduction is often utilized for 
the kind of systems studied in Sect. 5.2 and 
Sect. 5.3 (see, e.g., Antoulas [2005], Obinata & 
Anderson [2001]). 

5.4 An unstable monovariable nonlinear system 

5.4.1. i-PID For y — y = u'^ we utilize for Eq. 
(2) the local model y = F + u. The stabilization 
around a reference trajectory y* is provided by 
the i-PI controller 

u=-F + y* +Kpe + Ki j e (11) 

where Kp = —2, Kj = —1. The simulations 
displayed in Fig. 7 are excellent. 

5.4.2. Anti-windup We now assume that u 
should satisfy the following constraints — 2 < zt < 
0.4. The performances displayed by Fig. 8 are 
mediocre if an anti-windup is not added to the 
classic part of the i-PI controller. Our solution 
is elementary ^ : as soon as the control variable 
gets saturated, the integral / 6 in Eq. (11) is 
maintained constant. 



5.5 Ball and beam 

Fig. 10 displays the famous ball and beam exam- 
ple, which obeys to the equation y = Byii^ — 
BGsmu, where w = ^ is the control variable. This 
monovariable system, which is not linearizable by 
a static state feedback, is therefore not flat. It is 
thus difficult to handle. 

We have chosen for Eq. {2) y = F + lOOu. In 
order to satisfy as well as possible the experimen- 
tal conditions, the control variable is saturated: 
— 7r/3 < u < 7r/3 and — tt < ii < tt. Fig. 11 and 
Fig. 12 display two types of trajectories: a Bezier 
polynomial and a sine function. We obtain in both 
cases excellent trackings thanks to an i-PID con- 
troller. In the Figures ll-(b), 12-(b), ll-(c), 12- 
(c) the control variable and the estimations of F 



'^'^This is a well covered subject in the literature (see, e.g., 
Bohn & Atherton [1995], Hippe [2006], Peng, Vrancic & 
Ilanus [1996]. 

^® Better performances would be easily reached, as flatness- 
based control is teaching us, with a modified reference 
trajectory. 

^^The sine function which appears in that equation takes 
us outside of the theory sketched in Sect. 2. This difficulty 
may be easily circumvented by utilizing tg-| (see Fliess, 
Levine, Martin & Rouchon [1995[). 

A quite large literature has been devoted to this example 
(see, e.g., (Fantoni Lozano [2002[, Hauser, Sastry & 
Kokotovic [1992], Sastry [1999[) for advanced nonlinear 
techniques, and (Zhang, Jiang & Wang [2002]) for neural 
networks). Let us add that all numerical simulations in 
those references are given without any corrupting noise. 



are presented in the noiseless case. Compare with 
the Figures 11- (d) and 12-(d)) where a corrupting 
noise is added. 



5.6 The three tank example 

The three tank example in Fig. 13 is quite popular 
in diagnosis. It obeys to the equations: 

xi = -Cisign(a;i - X3)^\xi - Xa] +Ui/S 
±2 = C3sign(a;3 - X2)\/\x3 - X2\ 
-C2sign(a;2)-v/|a;2| + U2/S 
±3 = Cisign(a;i - X3)y^i ' 



X3\ 



-C3sign(a;3 - X2)y/\x3 - X2I 



2/1 



Xi 



2/2 = X2 
{y3= X3 

where 

Cn = {l/S).^in-Spy/2g,n= 1,2,3; 

S — 0.0154 m (tank section); 

Sp = 5.10~^ m (pipe section between the tanks); 

g = 9.81 m.s~^ (gravity); 

/^i = A*3 = 0.5, IJ.2 = 0.675 (viscosity coefficients). 

As often in industry we utilize a zero-hold control 
(see Fig. 14- (c)). A decoupled Eq. (7) is employed 
here as we already did in Sect. 5.3: i/i = Fi + 
200(i;, i = 1,2. Fig. 14- (a) displays the trajectories 
tracking. The derivatives estimation in Fig. 14- 
(b) is excellent in spite of the additive corrupting 
noise. The nominal controls (Fig. 14-(c)) are not 
very far from those we would have computed with 
a flatness-based viewpoint (see Fliess, Join & Sira- 
Ramirez [2005]). We also utilize the following i-PI 
controllers 



1 

200 



y* -Fi + lOei + 2.10- 



1,2 



where y* is the reference trajectory, = y*—yi. In 
order to get a good estimate of Cj we are denoising 
Ui (see Fig. 14- (d)) according to the techniques of 
Sect. 3. 



6. CONTROL WITH A RESTRICTED 
MODEL 

6.1 General principles 

6.1.1. Flatness The system (4) is assumed to 
be square, i.e., m = p, and flat. Moreover y = 



(2/I) ■ • ■ , 2/m) is assumed to be a flat output, i.e., 
Mj = 0, j = 1, . . . , m. It yields locally 

«j=fj(y,...,y(^^-)), i = l,...,m (12) 

Flatness-based control permits to select easily an 
efficient reference trajectory y* to which corre- 
sponds via Eq. (4) and Eq. (12) the open loop 
control u* . Let e = y — y* be the tracking error. 
We assume the existence of a feedback controller 
"feedback (e) such that 

u = u* + Mfeedback(e) (13) 

ensures a stable tracking around the reference 
trajectory. 



6.1.2. Intelligent controllers Replace Eq. (4) by 
$,(y,...,y(^^),M,...,u(*^-)) + G,=0 (14) 

where the Gj, j = 1, .... m, stand for the unmod- 
eled parts. Eq. (12) becomes then 



. ,m (15) 



where Hj ^ Gj in general. Thanks to Eq. (15), Hj 
is estimated in the same way as the fault variables 
and the unknown perturbations are in (Fliess, 
Join & Sira- Ramirez [2008]). Consider again y* 
and u* as they are defined above. The intelligent 
controller, or i-controller, follows from Eq. (13) 



u 



u* + 



(Hi 



"feedback (e) 



It ensures tracking stabilization around the refer- 
ence trajectory. 



6.2 Frictions and nonlinearities 

A point mass m at the end of a spring of length y 
obeys to the equation 



my = -K{y) + T{y) - dy + F^, 



(16) 



where 

• ^ext = u is the control variable; 

• d and F{y) are due to complex friction phe- 
nomena; 

• /C(y) = kiy + fcsj/^ exhibits a cubic nonlin- 
earity of Duffing type; 

The mass m = 0.5 is known; there is a possible 
error of 33% for fci = 3 and we utilize ki = 2; d 
and fcs, which are unknown, are equal to 5 and 10 
in the numerical simulations. For the frictions, 



See (Fliess, Join & Sira- Ramirez [2005]) for details and 
references. This paper was presenting apparently for the 
first time the diagnosis, the control and the fault accom- 
modation of a nonlinear system with uncertain parameters. 
See also Fliess, Join & Sira- Ramirez [2008]. 



There is a huge literature in tribology where various 
possible friction models are suggested (see, e.g., in control 
(Olsson, Astrom, Canudas de Wit, Gafvert Sz Lischinsky 
[1998], Nuninger, Perruquetti & Richard [2006]). Those 
modelings are bypassed here. 



we have chosen for the sake of computer simula- 
tions the well known model due to Tustin [1947]. 
Fig. 15-(a) exhibits its quite wild behavior when 
the sign of the speed is changing. 

6.2.1. A classic PID controller The PID con- 
troller is tuned only thanks to the restricted model 
my = —kiy + u. Its gains are determined in such 
a way that all the poles of the closed-loop system 
are equal to —3: Kp = —ki + 27m, Kj = —27m, 
Kd = 9m. 

6.2.2. The corresponding i-PID controller Pick 
up a reference trajectory y* . Set 

u* = my* + kiy* 

Our i-PID controller is given by 

Fext = w = w*- [a]e + PID(e) (17) 

where 

• Q = T{y) — {ki — ki)y — k^y^ — dy, which 
stands for the whole set of unknown effects, 
is estimated via 

[Q]e = m[y]e + k[y]e - Fe,t 

which follows from Eq. (16) {[y]e and [y]e are 
the denoised output variable and its denoised 
2"''-order derivative - see Fig. 15-(d,f)); 

• PID(e), e = y — y*, is the above classic PID 
controller. 

6.2.3. Numerical simulations The performances 
of our i-PID controller (17), which are displayed 
in Fig. 15-(c,d), are excellent. When compared to 
the Figures 

• 15-(e,f), where 

• flatness-based control is employed for 
determined the open-loop output and 
input variables, 

• the loop is closed via a classic PID con- 
troller, which does not take into account 
the unknown effects; 

• 15-(g,h), where only a classic PID controller 
is used, without any flatness-based control; 

the superiority of our control design is obvious. 
This superiority is increasing with the friction. 



6.3 Non-minimum phase systems 

Consider the transfer function 
s — a 
s'^ — {b + c)s + be 

a,b,c G M are respectively its zero and its two 
poles. The corresponding input-output system is 



non-minimum phase if a > 0. The controllable 
and observable state-variable representation 

Xl = X2 

< X2 = {b + c)x2 — bcxi + u (18) 

y = X2 — axi 

shows that z = cci is a flat output. The flat 
output is therefore not the measured output, as 
we assumed in Sect. 6.1.1. Our control design has 
therefore to be modifled. 

6.3.1. Control of the exact model To a nominal 
flat output z* corresponds a nominal control 
variable 

u* = z* - {b + c)z* + bcz* 

and a nominal output variable 

y* = z* - az* (19) 

Introduce the GPI controller (Fliess, Marquez, 
Delaleau & Sira- Ramirez [2002]) 

u= u* + -/J{u-u*)+Kp{y-y*) 

+KjJ{y-y*) + KuJJ{y-y*) ^^^^ 

where the coefScients j, Kp, Kj, Kn e M are 
chosen in order to stabilize the error dynamics 
e — z — z*. Excellent performances are displayed 
in Figures 16- (a) and (b), where o = 1, 6 = — 1, 
c = —0.5, even with an additive corrupting noise. 
See Figures 16-(c) and (d) for y* and z* which is 
calculated by integrating Eq. (19) back in time. 

6.3.2. Unmodeled effects The second hne of Eq. 
(18) may be written again as 

X2 = {b+ c)x2 — bcxi + u + m 

where zu stands for the unmodeled effects, hke 
frictions or an actuator's fault. Replace the nom- 
inal control variable u* of Sect. 6.3.1 by 

where [zuje is the estimated value of vu, which is 
given by 

He = - + Uj 

Moreover, 

[u]e= ii* +j{u-u*)+Kp{[y]^~y*) 
+Kji[yU-y*)+KiiJ{[yl-y*) 

follows from Eq. (20). 

Start with the Figures 17 and 18, where w = —0.5, 
a = 1, 6 = -1, c = -0.5. Fig. 17-(b), where the 
nominal control is left unmodifled, and Fig. 17- 
(e), where it is modifled, demonstrate a clear-cut 



^^This Section, which is based on previous studies (Fliess 
& Marquez [2000], Fliess, Marquez, Delaleau & Sira- 
Ramirez [2002]), should make the reading of Sect. 6.3.2 
easier. 



superiority of our approach, even with an additive 
corrupting noise. 

In Fig. 19, zu is no more assumed to be constant, 
but equal to — O.ly. In the numerical simulations, 
o = 2, 6 = — 1 c = 1. If n7 is not estimated, 
it influences the tracking quite a lot even if its 
amplitude is weak. When vu is estimated on the 
other hand, the results are excellent, even with an 
additive corrupting noise. 

7. CONCLUSION 

The results which were already obtained with 
our intelligent PID controllers lead us to the 
hope that they will greatly improve the practical 
applicability and the performances of the classic 
PIDs, at least for all finite-dimensional systems 
which are known to be non-minimum phase within 
their operating range: 

• the tuning of the gains of i-PIDs is straight- 
forward since 

• the unknown part is eliminated, 

• the control design boils down to a pure 
integrator of order 1 or 2; 

• the identification techniques for implement- 
ing classic PID regulators, which are often 
imprecise and difficult to handle, are becom- 
ing obsolete. 

Model-free control and the control with a re- 
stricted model seem to question the very princi- 
ples of modeling in applied sciences, at least when 
one wishes to control some concrete plant. This 
might be a fundamental "epistemological" change, 
which needs of course to be further discussed and 
analyzed. A natural extension to uncontrolled sys- 
tems is being developed via various questions in 
financial engineering: see already the preliminary 
studies in (Fliess & Join [2008b, 2009a,b]). 
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Figure 1. Stable linear monovariable system 
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Figure 2. Modified stable linear monovariable system 



(a) i-PI control 



(b) Output (-); reference (- -); denoised 
output (. .) 



(c) PID control 
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Figure 3. Stable linear monovariable system, with an actuator's fault 
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Figure 4. Linear monovariable system with a large spectrum 



(a) Output's noises 



(h) Estimation of y\ 



(c) Estimation of y2 





(a) i-PID: references (- -) and outputs (h) PID: references (- -) and outputs 

Figure 6. Linear multivariable system 
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Figure 7. Instable nonlinear monovariable system 
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Figure 8. Instable nonlinear system: saturated control without anti-windup 
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Figure 9. Instable nonlinear system: saturated control with anti-windup 




Figure 10. The ball and beam example 
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Figure 11. Polynomial trajectory for the ball and beam 




(a) Reference (- -); output 
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(c) Estimation of F (d) Reference (- -); output in the noisy case 

Figure 12. Sinusoidal trajectory for the ball and beam 
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ure 13. The 3 tank system 




( d) Denoising 

ure 14. Simulations for the 3 tank system 






(a) The friction law 



(h) Friction T(y) (-); Estimation of 
the whole unknown effects [5{y)]e 



(c) i-PID control 





(d) i-PID: output (-), reference (- -); (e) Flatness-based control 
denoised output (. .) PID 




classic (f) Flatness-based control + classic 
PID: output (-), reference (- -), de- 
noised output (. .) 





(g) Non-flatness-based control + PID (h) Non-flatness-based control: output 

(-), reference (- -), denoised output 
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Figure 15. The spring with unknown with nonhnearity, friction and damping 




Figure 16. Non-minimum phase system 



(a) u {-); u* (- -) 
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(d) u (-); u* (- -); Mpej.t(. .) (e) Output (-); reference (- -); denoised 

output (. .) 

Figure 17. The non-minimum phase system where the first effect is not modeled 




Figure 18. The non-minimum phase system where the first effect is not modeled 




Figure 19. The non-minimum phase system where the second effect is not modeled 
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